Poissonian reducibility and thermal scaling within the lattice gas model and molecular 

dynamics model 
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The emission of clusters in the nuclear disassembly is investigated within the framework of isospin 
dependent lattice gas model and classical molecular dynamics model. As observed in the recent 
experimental data, it is found that the emission of individual cluster is poissonian and thermal 
scaling is observed in the linear Arrhenius plots made from the average multiplicity of each cluster. 
The mass, isotope and charge dependent "emission barriers" are extracted from the slopes of the 
Arrhenius plots and their possible physical implications are investigated. 

PACS numbers: 25.70.Pq, 05.70.Jk, 24.10.Pa, 02.70.Ns 



In low-intermediate energy heavy ion collisions (HIC) 
the hot nuclei with moderate temperature can be formed 
and they finally deexcite by the different decay modes, 
such as the emission of multiple intermediate mass frag- 
ment (IMF), i.e. multifragmentation. Despite of the ex- 
tensive studies in experiments and theories, it is still dif- 
ficult to clarify whether the multifragmentation is statis- 
tical or dynamical, sequential or simultaneous. Recently 
Moretto et al. found that there exists the resilient re- 
ducibility and thermal scaling in multiple fragment emis- 
sion process, which gives one a helpful and clear picture 
to look and understand the multifragmentation. They 
observed that the experimental Z-integrated fragment 
multiplicity distributions P™ are binomially distributed, 
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in each transverse energy (E t ) window, where n is the 
number of emitted fragments and m is interpreted as the 
number of times the system tries to emit a fragment. The 
probability of emitting n fragments can be reduced to a 
single-particle emission probability p which gives linear 
Arrhenius plots (i.e. excitation functions) when ln(l/p) 
is plotted vs \/^/(E t ). By assuming a linear relationship 
between y/(Et) and temperature T, the linearity of the 
observed ln(l/p) vs l/y/(E t ) plot can be explained to 
a thermal scaling of the multifragment process [jl], |2|, || . 
In this case, these linear Arrhenius plots suggest that 
p has the Boltzman form p oc exp(—B/T) with a com- 
mon fragment barrier B. However, since the binomial 
decomposition has been performed on the Z-integrated 
multiplicities, typically associated with 3 < Z < 20, the 
Arrhenius plot generated with the resulting one fragment 
probability p is an average over a range of Z values. More 
recently, Beaulieu, Phair and Moretto et al. found that 
the fit with the binomial distribution can be also replaced 
with the Poisson distribution in the constraint of charge 
conservation. 

Instead of analyzing for Z-integrated multiplicities, 
they analyzed the behavior of individual fragment species 



of a given Z for higher resolution experimental data and 
noticed that the n-fragment multiplicities P(n) obey a 
nearly Poisson distribution, 



P(n) 



< n > 7 
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where n is the number of fragments of a given Z and the 
average value < n > is a function of the total transverse 
energy E t , and were thus reducible to a single- fragment 
probability proportional to the average value < n > for 
each Z Q. Similarly the < n > is found to be pro- 
portional to exp(—B/T) providing that T cx \jE t , i.e. 
there exists also a thermal scaling law. As pointed out by 
Moretto et al. ||, this kind of "reducibility" and "ther- 
mal scaling" are empirically pervasive features of nu- 
clear multifragmentation. "Reducibility" proves nearly 
stochastic emission process. "Thermal scaling" gives an 
indication of thermalization. More recently, Elliott and 
Moretto et al. discovered that the common features of 
Poissonian reducibility and thermal scaling can also be 
revealed in percolation and the Fisher droplet model || . 
Of course, we should keep in mind that the assump- 
tion of \/(E t ) oc T may be only valid for compound 
nuclei formed at low-to-moderate temperatures, but fail 
at higher temperatures in the experimental data M. In 
the present work, we will adopt the true temperature to 
study the reducibility and thermal scaling even though 
the assumption of y/(E t ) oc T is found to be also valid at 
low-to-moderate temperatures in the present model cal- 
culation. By investigating the mean cluster multiplicity 
as a function of temperature, we will illustrate that the 
Poissonian reducibility and thermal scaling is also valid 
for fragment emission in nuclear disassembly via the the- 
oretical reexamination of thermal equilibrium models. 

In this Letter, we will analyze the fragment multiplicity 
distributions for each individual fragment Z and A value 
in the framework of isospin dependent lattice gas model 
(I-LGM) and classical molecular dynamics (I-CMD). We 
will show that they are Poissonian and the associated 
mean multiplicities for each Z or A give linear Arrhenius 
plots as the experimental data illustrated in M. The A 
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and Z dependent barriers are extracted and investigated 
as a function of source size. Within our knowledge, this 
is the first time to explore the Poissonian reducibility 
and its thermal scaling for the individual fragment in 
the nuclear disassembly within the lattice gas model and 
molecular dynamics model. 

The lattice gas model was developed to describe the 
liquid-gas phase transition for atomic system by Lee and 
Yang ||. The same model has already been applied 
to nuclear physics for isospin symmetrical systems in 
the grand canonical ensemble Ipfl with a sampling of the 
canonical ensemble fn| [□], [lgtOI ||, P, [jj, 0, and 
also for isospin asymmetrical nuclear matter in the mean 
field approximation jt8). In addition, a classical molecu- 
lar dynamical model is used to compare with the results 
of lattice gas model. Here we will make a brief description 
for the models. 

In the lattice gas model, A (= N + Z) nucleons with 
an occupation number Si which is defined Si — 1 (-1) for 
a proton (neutron) or Si — for a vacancy, are placed 
on the L sites of lattice. Nucleons in the nearest neigh- 
boring sites have interaction with an energy e SiSj . The 
hamiltonian is written as 



A p2 

1=1 Kj 



(3) 



where Pi is the momentum of the nucleon and m is its 
mass. The interaction constant e SiS . is chosen to be 
isospin dependent and be fixed to reproduce the bind- 
ing energy of the nuclei Jig] : 



= O.MeV, 



c pi 1 



= -5.33MeV. 



(4) 



Three-dimension cubic lattice with L sites is used which 
results in pf = j^p of an assumed freeze-out density of 
disassembling system, in which po is the normal nuclear 
density. The disassembly of the system is to be calcu- 
lated at pf, beyond which nucleons are too far apart 
to interact. Nucleons are put into lattice by Monte 
Carlo Metropolis sampling. Once the nucleons have been 
placed we also ascribe to each of them a momentum by 
Monte Carlo samplings of Maxwell-Boltzmann distribu- 
tion. 

Once this is done the I-LGM immediately gives the 
cluster distribution using the rule that two nucleons are 
part of the same cluster if 



Pr/2p - £s zSj SiSj < 0, 



(5) 



where P r is the relative momentum of two nucleons and 
p is their reduced mass. This prescription is evidenced 
to be similar to the Coniglio-Klein's prescription |19| in 
condensed matter physics and be valid in I-LGM Jll|, |l2j, 
|l4| , |l6|| . To calculate clusters using I-CMD we propagate 
the particles from the initial configuration for a long time 
under the influence of the chosen force. The form of the 
force is chosen to be also isospin dependent in order to 



compare with the results of I-LGM. The potential for 
unlike nucleons is 



Vn P (r)(— < a) = A 
v np (r)(— > a) = 0. 
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In the above, ro = 1.842/m is the distance between the 
centers of two adjacent cubes. The parameters of the po- 
tentials are p = 2, q = 1, a = 1.3, B — 0.924, and A — 
1966 MeV. With these parameters the potential is min- 
imum at tq with the value -5.33 MeV, is zero when the 
nucleons are more than 1.3r apart and becomes stronger 
repulsive when r is significantly less than ro . The poten- 
tial for like nucleons is written as 



v pp (r)(r < r ) 
Vpp(r)(r > r ) 



v np (r) - v ap (r ), 

0. (7) 



This means there is a repulsive core which goes to zero 
at ro and is zero afterwards. It is consistent with the 
fact that we do not put two like nucleons in the same 
cube. The system evolves for a long time from the initial 
configuration obtained by the lattice gas model under the 
influence of the above potential. At asymptotic times the 
clusters are easily recognized. The cluster distribution 
and the quantities based on it in the two models can now 
be compared. In the case of proton-proton interactions, 
the Coulomb interaction can also be added separately 
and compared with the cases where the Coulomb effects 
are ignored. 

In this Letter we choose the medium size nuclei 129 Xe 
as a main example to analyze the behavior of individ- 
ual fragment emission during nuclear disassembly with 
the helps of I-LGM and I-CMD. In addition, the systems 
with A sys = 80 (Z sys = 33) and 274 (Z sys = 114) are 
also studied to investigate the possible source size depen- 
dence. In most case, pj is chosen to be about 0.38 po, 
since the experimental data can be best fitted by pf be- 
tween 0.3po and 0.4/?o in the previous LGM calculations 
O, E0| , which corresponds to 7 3 cubic lattice is used for 



Xe, 6 d for A s 



80 and 9 A for A s 



274 . In the 



condition of the fixed freeze-out density, the only input 
parameter of the models is the temperature T. In the 
I-LGM case, pf can be thought as the freeze-out den- 
sity but in the I-CMD case pf is, strictly speaking, not a 
freeze-out density but merely defines the starting point 
for time evolution. However since classical evolution of a 
many particle system is entirely deterministic, the initial- 
ization does have in it all the information of the asymp- 
totic cluster distribution, we will continue to call pf as 
the freeze-out density. 1000 events are simulated for each 
T which ensures enough statistics. 

One of the basic characters of the Poisson distribution 
Eq.(2) is the ratio a^J < ni >^ 1 where a^. is the vari- 
ance of the distribution and < rtj > is the mean multi- 
plicity. The first step we will show is this ratio. We give 
these ratios for clusters classified with different masses 
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FIG. 1: The ratio of Ui 2 / < rii > for the clusters classified with 
mass, light isotope mass and atomic number as a function of tem- 
perature. The left panel is for the l-LGM calculation and the right 
for I-CMD with Coulomb. The symbols are illustrated on the fig- 
ure. 

(A), light isotopes (ISO) and atomic numbers (Z) for 
the disassembly of 129 Xe as a function of temperature in 
the framework of l-LGM and I-CMD with Coulomb in 
Figure 1. Obviously, the ratios are close to one except 
for protons, which indicates that they might belong to 
the Poisson distributions. This can be further supported 
by the multiplicity distribution in different temperature. 
For instance, Fig. 2 show the quality of the Poisson fits to 
the charged particle multiplicity distribution for 129 Xe 
in the l-LGM case. These Poisson fits are excellent for 
all Z > 2 over the entire range of T. The same good 
Poisson fit is obtained in the cases of I-CMD. Thus we 
can conclude that Poissonian reducibility is valid in the 
thermal-equilibrium lattice gas model or molecular dy- 
namics. 



To verify thermal scaling in the models, the temper- 
ature dependence of the mean yield of clusters is inves- 
tigated. Consequently, we generate Arrhenius plots by 
plotting In < nz > vs 1/T. Figure 3 give a family of 
these plots for the disassembly of 129 Xe within the frame- 
work of l-LGM (left panel) and I-CMD with Coulomb 



FIG. 2: The excitation functions P n for elements with Z > 2 
emission from the source 129 Xe in the l-LGM calculation. The 
lines are Poisson fits with Eq. (2). 



interaction (right panel). The observed Arrhenius plots 
are striking linear for the lower T side, and their slopes 
generally increase with increasing Z value. Generally, 
the thermal scaling is expected when the yields, for a 
fixed nucleon number system, are dominated by fragment 
binding. This is the case when the temperatures are low 
compared to the binding energy per particle. At these 
temperatures, one can anticipate one large fragment sur- 
rounded by many small clusters. The contrary tendency 
reveals in the high T side where In < nz > increases 
with 1/T, i.e. decreases with increasing T. In this case, 
nuclear Arrhenius plots of < n > with 1/T are not valid 
but the Poissonian reducibility still remains (see Fig. 2). 
This behavior of < n > at higher T is related to the 
branch of the fall of the multiplicity of IMF (Nimf) 
with T where the disassembling system is in vaporiza- 
tion [^2| ^||, HJj and hence only the lightest clusters are 
dominated and the heavier clusters become fewer and 
fewer with increasing T . Afterwards we will focus on the 
branch of lower temperature to discuss the Arrhenius law. 
The overall linear trend illustrates that thermal scaling is 
also present when the individual fragments of a specific 
A, ISO and Z are considered. 



4 



From figure 3 the slope parameter can be directly ex- 
tracted in the lower T side as a function of Z or A. In Ref. 
j| Moretto et al. has interpreted these slope parameters 
as "emission barriers" of specific individual fragments. 
Figure 4 gives the emission barrier of individual frag- 
ments with different A, ISO and Z in the framework of 
I-LGM, I-CMD with/without Coulomb interaction. The 
error bar in the figure represents the error in the extrac- 
tion of the slope parameter. The first indication from this 
figure is that the emission barrier in the FLGM case is the 
nearly same as the FCMD case without Coulomb force, 
which supports that I-LGM is equivalent to I-CMD with- 
out Coulomb interaction rather well when the nuclear 
potential parameter is moderately chosen, but I-LGM is 
a quick model to analyze the behavior of nuclear dissoci- 
ations. The inclusion of long-range Coulomb interaction 
makes the emission barrier of individual fragments much 
lower since the repulsion of Coulomb force reduces the 
attractive role of potential and hence make clusters es- 
cape easily. The second indication is that the emission 
barriers increase with A (Z) at low A (Z) values and tend 
to be saturated at high A (Z) ones. Similar experimental 
results have been observed for individual fragments with 
different Z in Ref. Q or different A in Ref.|(| . However, 
the middle panel of Fig. 4 shows that bare dependence 
of emission barrier of ISO on A in the fixed atomic num- 
ber Z, which indicates that the Z dependence of barrier 
is perhaps more intrinsical the A dependence is mostly 
due to the average effect over the species with the same 
A but different Z. 

I 



for the particles classified with different mass number, or 



for the particles classified with different charge number, 
where C2 is a fit constant for surface term and rcoui is 
chosen to be 1.22 fm. Z res is a fitted average charge 
number of the residue. (A/Z)fi t is chosen to be 1.866, as 
taken from the fits for I-LGM. The overall fits for A and 
Z dependent barrier in the case of I-CMD with Coulomb 
force give C2 = 12.921 and Z res ~ 41 with the dot-dashed 
line in Fig. 4a and 4c. The excellent fit supports that 
the Coulomb energy plays another important role in the 



On the origin of these barrier, the surface energy and 
Coulomb energy would play the roles. If the cluster emis- 
sion is mainly controlled by its surface energy, it would 
suggest barriers proportional to Z 2 / 3 (A 2 / 3 ). In the case 
of I-LGM and I-CMD without Coulomb, we can try to fit 
the barrier for the particles with different mass number 
by 

5coui.o// = cixA, 2/3 , (8) 

or for the particles with different charge number by 

Bcoui.off = ci x {(A/Z) flt * Zif/ 3 , (9) 

where (A/Z) fu is a fit coefficient of A/Z for emitted par- 
ticles, and Ai (Zi) is the mass (charge) of particle. c\ is 
the fit constant for surface energy term. The solid line 
in the Fig. 4a is a function of Eq.(8) with c\ = 8.469 and 
the solid line in the Fig. 4c is a function of Eq. (9) with c\ 
= 8.469 and (A/Z) fit = 1.866. These excellent fits imply 
that the surface energy play a major role in controlling 
the cluster emission when the long range Coulomb force 
is not considered. However for the cluster emission with 
the Coulomb field, we can assumed that the barrier is 
mainly constituted by the surface energy term and an 
additional Coulomb term as 



(10) 



(11) 

I 

cluster emission. 



In the case of I-LGM and I-CMD without Coulomb, 
one would expect the barrier for each Z (A) to be nearly 
independent of the system studied if only the surface en- 
ergy is substantial to the emission barrier. The left panel 
of the figure 5 shows the results for B^, Bjgo arj d Bz 
for three different systems in the I-LGM case. The same 
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FIG. 3: The average yield per event of different clusters classified 
with A (top), ISO (middle) and Z (bottom) as a function of 1/T. 
The left panel is for the I-LGM calculation and the right for I-CMD 
with Coulomb. The solid lines are fits to the calculations using a 
Boltzmann factor for < Hi >. The symbols are illustrated on the 
figure. 



FIG. 4: The emission barrier extracted from the Arrhcnius plots 
as a function of cluster mass (top), isotopic mass (middle) or clus- 
ter charge (bottom) in the cases of I-LGM (solid squares), I-CMD 
without Coulomb (solid circles) and with Coulomb (open circles). 
The solid lines are fits with the Eq. (8) or (9), and the dot-dashed 
lines represent the fits with the Eq. (10) or (11). 



freeze-out density of 0.38po and the same N/Z is chosen 
for the systems of A sys = 80 and A sys = 274. Actually, 
it appears to have no obvious dependence of emission 
barrier on source size as expected for the role of surface 
energy. The solid line in the figure is the same as in 
Fig. 4. However, when the long-range Coulomb interac- 
tion is considered, the emission barrier reveals a source 
size dependence. The right panel of figure 5 gives the 
emission barriers Ba 1 Bjso and Bz in the case of I-CMD 
with Coulomb force. It looks that the barrier increase 
with the decreasing of charge of system, which can be 
explained with the Eq. (10) and (11) where the decreas- 
ing of the residue Z res will result in the decreasing of the 
Coulomb barrier and hence the increasing of the emission 
barrier. The lines represent the fits with the Eq.(10) and 
(11) for three different mass systems. 



In the above calculations, the freeze-out density of sys- 
tems is fixed at <~ 0.38po- Considering the freeze-out den- 
sity is an important debating variable in the latter stage 
of heavy ion collisions, here we will discuss the possi- 



ble influence of freeze-out density on the emission barrier 
of clusters. The calculations at the freeze-out density 
of 0.177po and 0.597po f° r 129 Xe, corresponding to 9 3 
and 6 3 cubic lattices respectively, are supplemented to 
compare. Figure 6 gives the results of Ba, Biso and 
Bz at different density. It looks that there are no ob- 
vious freeze-out density dependence in the both cases of 
I-LGM and I-CMD. This is also consistent with that as- 
sumption that the surface energy is the dominant role in 
controlling the cluster emission. 



In conclusion, the poisson reducibility and thermal 
scaling of the emitted clusters is explored in the lattice 
gas model and molecular dynamical model. The calcula- 
tions are qualitatively consistent with the recent exper- 
imental observation by Moretto/Wozniak's group even 
though the temperature is supposed to be proportional 
to the total transverse energy in the latter experiments. 
A systematic study of the emission barrier on the cluster 
mass, isotope and charge proves that the cluster emis- 
sion is mainly controlled by both the surface energy and 
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FIG. 5: The source size dependence of the emission barriers for the 
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or cluster charge (bottom) from in the cases of I-LGM (left panel), 
I-CMD with Coulomb (right panel). The lines in the left panel are 
fits with the Eq. (8) or (9), and the solid, dot-dashed and dotted 
line in the right panel represents the fits to the emission barrier of 
A sys = 80, 129 and 274, respectively, with the Eq. (10) or (11). 
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the Coulomb interaction. In the framework of the lat- 
tice gas model and molecular dynamics model without 
the Coulomb interaction, the emission barrier relies on 
the cluster charge with the Z 2 / 3 (A 2 / 3 ) law and it docs 
not depend on the the source size and freeze-out den- 
sity which indicates that the surface energy play a basic 
dominant role to control the cluster emission. Conversely, 
in the framework of molecular dynamics model with the 
Coulomb force, the emission barrier will decrease strongly 
according to the Eq.(10) and (11) and it decreases with 
the increasing of the source size, illustrating that the 
Coulomb interaction also play another weighty role to 
control the cluster emission. 



FIG. 6: The emission barrier of 129 Xe for the different clusters 
classified with mass (top), isotopic mass (middle) or cluster charge 
(bottom) at the different freeze-out density in the cases of I-LGM 
(left panel), I-CMD with Coulomb (right panel). The lines are fits 
with the Eq. (8) or (9) in the left panel and with the Eq. (10) or 
(11 ) in the right panel. 
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